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ABSTRACT 


A recent innovation for determining the set of governing differen- 
tial equations of motion of a complex structure is described. Numerical 
values for the mass, stiffness, and damping coefficients of the dynamical 
equations associated with a particular input— response or transmission 
path are computed from the data usually obtained in conventional vibra- 
tion tests of a structure. The theory is based on the dynamic properties 
of multi-degree-of— freedom linear systems. The method requires the 
steady response (acceleration, velocity, displacement, or stress) and 
the driving sinusoidal force input for transmission paths of interest to 
be experimentally determined for a few frequencies near each major struc- 
ture resonance. Application of the method is illustrated by determining 
from experimental data the equations of motion of l/lO— and l/40— scale 
models of the Apollo /Saturn V launch vehicle. Transient responses computed 
using the identified equations for the 1 /40-scale model are shown to agree 
favorably with experimental results. 


INTRODUCTION 


The purpose of this paper is to present a simple technique for experi- 
mentally determining an acceptable set of equations of motion for a space 
vehicle structure. This system identification approach, while yielding 
the equations of motion, does not involve a detailed analysis or physical 
idealization of the structure. It relies solely on the experimental deter- 
mination of the steady state response of the structure to a sinusoidally 
varying input force for a few frequencies near each important structural 
resonance and on the usual assumptions regarding the behavior of lightly 
damped, linear structures. 



SYMBOLS 


Cm] 

tc] 

[K] 

ip(t)3 

Cx} 

*k 

*k 

£x}<J) 

x(J) 

k 

m 

icp} ( «J ) 

<n( <3 ) 


system mass matrix} symmetric, positive definite 
system viscous damping matrix} symmetric, nonnegative definite 
system stiffness matrix} symmetric, nonnegative definite 
system forcing function vector 
system coordinate vector 

system coordinate} an element of lx} 
amplitude of steady state displacement response at k 
system coordinate vector due only to response in the jth mode 
kth element of 

matrix of modal vectors {cp}( j) 

as columns 

modal vector obtained from |[K] - Xj[M^{x} = 0 
the ith and kth elements of {9} 


[*CJ 

C*KJ 


Cq} 

q o 

Fi(t) 

F. 


m(J),c ( ^ 
ik ' ik 


ik 

(J) 

ik 


a) 


diagonal modal mass matrix 

diagonal modal damping matrix 

diagonal modal stiffness matrix 

normal coordinates defined by tx} = [¥]tq} 

jth element of [q] 

jth modal mass, damping and stiffness 

arbitrary external point force applied at i 

amplitude of F-^siniut 

( 

,k^" / effective mass, damping and stiffness for the i-k 
input-output path in the jth mode 

phase angle between response and input force for the i-k path 

effective percent of critical damping in the jth mode 

circular frequency 


2 



UDj jth resonant circular frequency 

Units the units used in this paper are lb, in, sec 

THEORY 

The starting point for the development of the identification technique 
of this paper is in the mechanics of lightly damped, linear systems. The 
class of structures for which the technique is proposed is assumed to possess 
the following somewhat qualitative features which are stated in the form of 
assumptions and are basic to the ensuing development: 

1. Light damping typical of © space vehicle with no damping speci- 
fically designed into the structure. 

2. The modes of interest are sufficiently uncoupled in the velocity 
terms and separated in frequency so that a single degree-of— freedom 
analysis is adequate to represent the steady response above the 
half power point in a mode of interest. 

These basic assumptions are designed to imply that the steady response 
in each vibration mode of interest is not significantly affected by any 
other mode and that each mode can be isolated and individually exploited as 
discussed in detail in references (l) — (4). 

With the above assumptions as a background, the proposed identification 
technique and the supporting analytical arguments are now developed. 

General . — The equations of motion for the system may be written (ref. 5) 

[ M]&} + [C]&} + [K]{X} = tP(t) } (1) 
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The transformation 


( 2 ) 


fe} - min) 

■where 

\ (5 > 

is assumed to produce uncoupled equations in terms of the normal coordi- 
nates, q, so that 

nyai + csjtt} + ricjtq} - m T iFwi « 

where 

c"o,3 - ['n T cc]m 
c"k, 3 = m T wm 


The equation for the jth normal coordinate is 


+ c 0^j + k i q j “ 2cpp)F i (t) 


( 5 ) 


For only one external forcing function applied at i equation (¥) becomes 

( 6 ) 


m 


■J«J + + * ‘4'’ )p i< t > 


Near gesonant response . - Equation (2) may be written 


£x} - S {tp}0)q 
J 3 


so that the response in the jth mode is 

lx}^ = 


( 7 ) 
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or 


*£ J) = <£ 3 S 


( 8 ) 


and 


9 J ‘Ji) 

k 


( 9 ) 


Substituting equation ( 9 ) into equation (6) and dividing by cp' 


ID 


*1 \-(J) 

°,1 \ 

•U) 

k 3 1 

^3J^5 Tj*k + 

^Cp( J) 

*k + 



4 3) - *.(*) do) 


or 


m ji )x k ) > + c ik^*k^ + ^ik^^ ■ ^iC*) 


(11) 


The significance of equation (11) is that it represents the system 
dynamics for response at point k due to forcing at point i. Also, the 
derivation of equation (11) has not placed any restrictions on the forcing 
function at i. F i (t) is, in fact, an arbitrary, external, point forcing 
function. When F^(t) = F^sinaot with u) = ud^ only one equation of the 
form of equation (11) is required to represent the system response. It is 
upon this fact that the identification scheme is based. 

Equating the coefficients of equations (10 ) and (11), it is evident 


that 


m 


3 ‘ m ii *i J) 4 3) 


Cj * <4«5) cpj^) 

k <5 “ 4k 4 3) 4 J) 


(12) 
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Setting i = k 




( 13 ) 


from -which the modal amplitudes at i and k are found to he given by 



m 5 

m(«)) 

ii 

M) - . 



( 1*0 


m' 


ik i 

The quantities mu and m^ are defined to be positive and hence <p(«)) 

is positive. A negative value for cp^ implies a negative m J which 

k ^ 

is wholly consistent with the facts. 

Solutions for m-y^ , cy^, and k^ k . — The steady state or particular 
solution of equation (11) with 'F^(t) * F-^sinurt is given by 

x - X sin(<ut ~ 0) (15) 

Substituting equation (15) into equation (11) and solving for the 
coefficients Cy^ and kyj results in 

Fi 


and 


k ik } = 4*^ + ^ cos6 ik 


c ik^ “ ~ sin9 ik 
ik uK k ** 


( 16 ) 


(17) 


If one set of values of iw, F i# X k and 9 ik for equation (17) and two 
sets for equation (16) are known, the coefficients of equation (11) may be 
computed. 
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Response to Arbitrary Force . — The total system response at k due 
to any time dependent force Fj.(t) at i is found by superposition of 
the solutions of each of the set of the following equations: 


+ + k ik ) 4 1) * F i (t) 




and is given using equations (3) and (8) as 


= F i (t) 


“ik^ + C ik } 4 j) 

• ! 

m (p) x (p)+ c^F)x^ p ^ + k^^x(.P) = fL ( t) 

ik k ik k ik He i v ' 


x,, = 2 


(J) 




j=l 


(18) 
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EXPERIMENT/!, PROCEDURE 


The Langley l/lO— scale and l/4o— scale models of the Apollo /Saturn V 
launch vehicle are shown in figures 1 and 2 and the coordinate systems for 
both models are presented in figure 3* The l/lO-scale model is fully 
described in reference 6. 

l/lO-scale model . — The l/lO-scale model was complete in the lift-off 
structural configuration, but was entirely empty of simulated propellants. 

The boundary conditions were cantilevered— free. A steady frequency, trans- 
verse, sinusoidally varying input force was applied in the pitch plane 
through a strain gage type force gage at station $86 and the displacements 
at stations 4l8 (the tip of the escape tower), 377, and 282 were measured 
using a contacting, cantilever, strain gage beam. The strain gage displace- 
ment transducers had natural frequencies in the order of 60 cps and were used 
within their flat response regime. The signals from the force and the 
displacement transducers were processed through a balancing bridge, diffe- 
rential amplifier, D.C. isolation amplifier and were then recorded together 
with the calibration signals on an FM analog tape recorder. The selection 
of both strain gage force and displacement transducers with both signals 
processed through identical electronics was to assure the accurate determina- 
tion of phase angle, 0 ik , for use in equations (l6) and (17) . 

lAo— scale model . — The lAo-scale model was also complete in the lift- 
off structural configuration. Propellant loading corresponded to first 
stage burnout and the boundary conditions simulated were free— free. A steady 
frequency transverse, sinusoidally varying input force was applied in the 
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pitch, plane through a crystal type force gage alternately at stations 0 
and k-2 and the acceleration responses were measured at model stations 
102.9 and -2.7. The crystal transducers were used with this model in 
order to evaluate the quality of data produced by the two different 
instrumentation schemes. The signals from the force gage and accelero- 
meters were processed through similar conditioning equipment to minimize 
relative phase shift and, together with calibration signals, were recorded 
on tape. 

Data reduction . — The experimental analog data were digitally filtered 
using a 2k — point per cycle Fourier analysis from which the numerical ampli- 
tudes of the fundamental components of the input force and of the displace- 
ment and acceleration responses and the input— response phase angles, 9 ik> 
were computed. 


IDENTIFICATION PROCEDURE 


System Equations 


For convenience equations (18) were written as 


..(1) o (1) .(1) • 2 (1) 

*k + 2 ^ik Vk +<U L X k 

I I I 





**(P) 
*k + 



! 

_ 1 


m’ 


(a p i (t) 

ik 


and 


*k 


P 

£ 

0=1 



(19) 


( 20 ) 
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Computer experiments , — A controlled computer experiment was 
conducted in which an exact numerical solution of equations (19) 
and (20) was generated for an assumed set of typical coefficients. 

The solution was then corrupted and the coefficients m^jp , c^ , 
and kjj'p were computed using equations (16) and ( 17 ). The results 
indicated that errors of up to + 10 percent in F^/x^ and + 10 
degrees in 9^ can he tolerated hy this identifier when applied to 
systems which satisfy the basic assumptions of this paper. Since it 
was felt that the accuracy of the experimental data fell within these 
limits, identification of the l/lO-scale and lAo— scale models was 
attempted. 

Solution for coefficients . — The digitized experimental data, F^, 

«• 

x k or , 9 ik corresponding to a near resonant value of «u were 
used to determine the values of m^) , cj^ , kjjjp for each signifi- 
cant mode for each of the selected input— response paths for each model. 
Four or five frequencies near each resonance were used. The value of 
C (J) for each point was computed using equation (17)* The values of 

IK 

m(«5) and k(3) were determined hy solving equation (16) as a pair of 
simultaneous equations for slightly different near-resonant values of 
w. The average value of effective mass, damping, and stiffness for each 
mode was then determined. Typical first mode values for station 
577 response of the l/lO-scale model are as follows: 
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Frequencies 

Coefficients of Eqn. 

Coefficients of Eqn. 


f 

U) 




ik 

•5-sft 

ik 


4.52 

28.39 

• 34-79 

.1458 

304.9 

.4191 

876.4 

2.87 

4.72 

29.64 

• 3381 

.1615 

296.2 

.4777 

876.7 

2.96 


31-34 

• 3311 

.1611 

290.I 

.4866 

876.2 

3.02 


The sign of m^) computed using equation (.16) was very simply verified. 


ik 


The sign was taken to be positive if was in phase (approximately 0°) 

with cpjp) and negative if cp( was out of phase (approximately l8o°) 

with cp( «3) » The identification results for both models are listed as coeffi— 
i 

cients of equation (19) in Tables I and II. 


RESULTS AND DISCUSSION 


Comparisons of the identification and experimental frequency responses 
for the two models are presented in figures 4—10. 

l/lO-scale results.- The results for the l/lO-scale model are shown in 


figures 4-6, plotted as the response ratio 


Xr 


versus frequency. The 


particular solutions of equations (19) were obtained for F i (t) = sinuut, 

i = station 386, for 0 £ — ^30 using the coefficients listed in Table I. 

2n 
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The frequency response at stations 4l8, 377, and 282 were then obtained 
using equation (20)« For example, 

*418 ' *£!b + *418 + 4ll 

lAo-scale results . — The results for the l/4o— scale model are shown 
in figures 7—10. The frequency response solutions for forcing first at 
station 0 and then at station 42 were computed using equations (19) and 
(20) and the coefficients given in Table II, For example, for response at 
station 102.9 due to forcing at station 42 

*102.9 + 5 - l4i io2.9 + 69 ’ 530 4JL 9 = - 10 ’ 1 Si "“* 

*102. 9 + 4 - 54 *lo2. 9 + 289 ’ 300 *1 o 2. 9 ’ l6 ' 2 SimUt 

x (5) + 13.71x^ 5 ^ + 68l,l6ox^ 5 ! = -12.2 sintut 

102.9 102.9 102.9 

•<4) „ .(4) „ (4) 

X 102.9 + 26 * 25X 102.9 + 929,8T0X 102. 9 “ ~ n * 7 Sin<wt 


and 


x - x* 1 ) + x( 2 ) + x(3) + x (^) 

102.9 102.9 102.9 102.9 102.9 


Discussion. — The identification results for both models agree quite 


favorably with the experimental response data. In addition, the associated 


phase angles. 



, also agree, usually to within 5 to 10 degrees. 
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Therefore, it is felt that, for both models and for the input— response 
paths investigated, systems of equations suitable for computing the 
response to an arbitrary forcing function have been obtained. 

In this connection, the uniqueness of the identified equations has 
not been rigorously established. It is, however, felt that amplitude 
agreement as shown in figures 4—10 together with phase agreement to within 
5 to 10 degrees, constitutes sufficient conditions for a system identifi- 
cation adequate for all engineering purposes. The experience thus far 
also indicates, for the class of systems considered in this paper, that 
an identifier based on requiring coincidence of frequency response ampli- 
tudes only, without regard to phase, will produce basically the same 
results as if phase information were employed. This fact can be useful, 
because, in general, some obvious small adjustments of the identified 
parameters to give better results based on amplitude comparison is usually 
possible. Figures 4—10 as presented indicate the results that were achieved 
using the method of this paper without iteration. However, parameters were 
easily selected to produce perfect coincidence of the amplitude plots of 
both figures 8 and 10, for example. The phase differences were negligible. 
Therefore, it is felt that refinement of the identified equations to pro- 
duce perfect amplitude agreement is permissable, if not desirable. 

The advantages of using the approach of this paper are that a detailed 
structural idealization and associated analytical model is not required. 

The modes that actually contributed to the response were immediately iden- 
tified as the only ones observable for a given input response path. This 
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obviated the usual concern over the problem of including all of the 
significant vibratory modes of the structures. 

Once the equations of motion for a structure have been identified 
the transient response to an arbitrary force can be computed with confi- 
dence. In this connection, rigid body modes, if required, can be 
calculated from model drawings or experimentally determined. Results 
thus far of transient tests with the l/4o~scale model have produced 
excellent agreement between the acceleration response computed using the 
identified equations and the experimental transient acceleration responses. 
For example, a comparison of identification and experimental transient 
acceleration and displacement responses for the l/4o-seale model is shown 
in figure 11(a) and (b). The coefficients used in computing the accelera- 
tion response as predicted using the system identification results are 
given in Table II for i *» 0 and k = 102.9 (see figure 7)® In addition, 
the two rigid body modes required for this free— free system were included 
which resulted in five equations of motion which were integrated using the 
measured values of the transient input shown in figure 11(c) which was 
applied at station 0 of the model. The identified equations predict 
transient acceleration and displacement responses at station 102.9 which 
are in good agreement with the experimentally determined response. The 
results of figure 11 are especially important because they demonstrate the 
ability of the identified equations to predict accurately the response to 
an input of a different character than was used for their derivation. 
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CONCLUDING REMARKS 


A technique for determining the equations of motion of a complex 
structure has been presented. Both the number of the essential degrees 
of freedom and the coefficients of the equations are determined by the 
procedure which is applicable to a large class of aerospace and other 
structures. The procedure requires that good quality experimental frequen- 
cy response data be obtained for the significant resonances associated 
with specified input response paths. 

It is felt, for the class of structures considered in this paper, that 
coincidence of the identified and experimental frequency response amplitudes 
constitutes a sufficient condition for a satisfactory identification. This 
hypothesis, if true, suggests that refinement of the initial identification 
results to produce perfect amplitude agreement is desirable and may be 
useful in the formulation of an identifier that does not require explicit 
experimental phase information. 

It is felt that this identification procedure should be advantageous 
when the dynamical equations of motion for an existing structure are desired. 
For example, if the available resources are not sufficient to allow the for- 
mulation and verification of a detailed analytical model of an existing 
structure, of if only sinusoidal test equipment is available but transient 
response data are required, the procedure of this paper may be very useful. 
In particular, the identified equations were shown adequately to predict the 
system acceleration response to an arbitrary transient force. 
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TABLE I.- IDENTIFIED VALUES OF THE COEFFICIENTS OF EQUATION ( 19 ) 

FOR l/lO-SCALE MODEL 
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TABLE II. - IDENTIFIED VALUES OF THE COEFFICIENTS OF EQUATION (19) FOE l/bO-SCMJE MODEL 
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Figure 2.- l/4o-scale Apollo/Saturn V model. 
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Figure 3.- Schematic of Apollo /Saturn V models. 



POINT OF RESPONSE k: STA. 418 



Figure 4 .- Comparison of identification results with experimental frequency response for 

l/lO-scale model. 
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Figure 6 • - Comparison of identification results with experimental frequency response for 

l/lO-scale model. 



POINT OF RESPONSE k: STA. 102.9 
INPUT FORCE i: STA. 0 



Figure 7*- Comparison of identification results -with experimental frequency response for 

l/40- scale model. 
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POINT OF RESPONSE k: STA. 102.9 
INPUT FORCE i: STA. 42 



Figure 9-- Comparison of identification results with experimental frequency response for 

l/40-scale model. 



POINT OF RESPONSE k : STA. -2.7 
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Figure 10.- Comparison of identification results with experimental frequency response for 

l/40-scale model. 



(a) ACCELERATION RESPONSE AT STATION 102.9 



NASA-Langley, 1968 


Figure 11.- Comparison of identification and experimental transient response for 

l/40- scale model. 



